#ifndef LINEAR_EQ_BAND_SYMM_MATRIX_H_
#define LINEAR_EQ_BAND_SYMM_MATRIX_H_

void solve_lower_system(double* matrix, int order, int half_bandwith, double* b);

void solve_upper_system(double* matrix, int order, int half_bandwith, double* b);

void solve_system_factorized_matrix(double* matrix, int order, int half_bandwith, double* b);


#endif /* LINEAR_EQ_BAND_SYMM_MATRIX_H_ */